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Abstract. We provide a brief characterization of the main features of the homo- 
geneous sine- Gordon models and discuss a general construction principle for colour 
valued S-matrices, associated to a pair of simply laced Lie algebras, which contain 
(N ' the homogeneous sine-Gordon models as a subclass. We give a brief introduction to 

the thermodynamic Bethe ansatz and the form factor approach and discuss explicit 
solutions for both methods related to the homogeneous Sine-Gordon models and its 
generalization. 
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1. Introduction 
y—i [ 

The work I want to report about |l], ^, ||, [|, |(| is based on a 
' ■^ ' collaboration with Olalla Castro-Alvaredo (Universidad de Santiago 

de Compostela), Christian KorfF (Stony Brook) and Luis Miramontes 
D I (Universidad de Santiago de Compostela). 

The completion of the entire bootstrap program [|?J for integrable 
quantum field theories in 1+1 space-time dimensions remains an open 
challenge for most concrete models. Roughly speaking "the program" 
can be divided into five distinct steps: i) the determination of the exact 
on-shell S-matrix, ii) the computation of closed formulae for the n- 
particle form factors, iii) the identification of the entire local operator 
content, iv) the computation of the related correlation functions and v) 
various non-perturbative consistency checks, like the thermodynamic 
Bethe ansatz, which exploit the information provided by the related 
conformal field theories. Step v) may of course also be carried out by 
resorting to standard old fashioned perturbative computations for vari- 
ous quantities. However, this is slightly opposed to the main virtue and 
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appeal of this approach, which is its entirely non-perturbative nature 
in the coupling constant. 

Recently we investigated |l], ||, ||, HI [5], |j a class of models, the homo- 
geneous Sine-Gordon model (HSG) |8| and its generalization, for which 
the bootstrap program was completed to a large extend. In comparison 
with other models, describable by a well-defined Lagrangian and which 
have been studied so far in 1+1 dimensions, the HSG-models possess 
two interesting novel features familiar from realistic 1+3 dimensional 
theories. Namely, their spectrum contains unstable particles and parity 
invariance is broken. With regard to one of the motivations for the 
study of lower dimensional theories, namely to treat them as test labo- 
ratories for general concepts and ideas, these properties deserve further 
consideration. One should say that unstable particles have first been 
investigated in 1+1 dimensions in the context of the so-called roaming 
(staircase) models @ and parity breaking was previously observed for 
the Federbush models in which the scattering matrix is however 
just a trivial constant, i.e. rapidity independent, phase. 

Due to the limited amount of time and space I will focus in my 
account on the main findings, referring the reader to the original lit- 
erature [|], ||, [3|, H], H ^] and |ll], 12, 13] for more details and general 
background. 



2. The homogeneous sine-Gordon models 

I want to start by recalling the main features of the G^-HSG-models || 
related to simply laced Lie algebras. The latter have been constructed 
as integrable perturbations of gauged WZNW-coset theories ]l4| of the 
form Gk/U(lY , where G is a compact simple Lie group of rank t and 
k > 1 the level of the Kac-Moody algebra. The defining action of the 
HSG-models reads 

ShsgM = SwzNwb] + — ^ / < ^ x ( A +>9(x)' A-g(x)) • (1) 

vrp J 

Here Swznw denotes the coset action, ( , ) the Killing form of G and 
g(x) a group valued bosonic scalar field. A± are arbitrary semi-simple 
elements of the Cartan subalgebra associated with the maximal abelian 
torus H C G, which have to be chosen not orthogonal to any root of 
G. The parameters m and (3 2 = 1/k + 0(l/k 2 ) are the bare mass scale 
and the coupling constant, respectively. 

The study of this theory from a classical point of view has been 
carried out already to a very large extend and constitutes meanwhile 
a subject in its own right. Here I just want to focus on the properties 
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which makes it an interesting theory from the quantum field theoretical 
point of view. In fact, as will be argued in the next sections, from that 
standpoint, one could abandon the classical picture altogether, since it 
just provides complementary but not necessary information. This is of 
course the natural algebraic quantum field theorist's viewpoint. 

The classical equations of motion of these models correspond to 
non-abelian affine Toda equations M, 15 1, which are known to be clas- 



sically integrable and admit soliton solutions. The integr ability on the 
quantum level was established in |l(| by the construction of non-trivial 
conserved charges. 

The coset conformal field theory defined by SwzNwb] is character- 
ized by data which are in principle extractable from the thermodynamic 
Bethe ansatz (TBA) and the form factor analysis, namely the Virasoro 
central charge c of the coset and the conformal dimensions A of various 



primary fields [14] 



a / a u (A-(A + 2p)) (A -A) 

c ^m e = kTh M > A(A ' A)= 2(k + h) 2T- (2) 

Here h is the Coxeter number, p the Weyl vector, i.e. the sum over 
all fundamental weights, A a highest dominant weight of level smaller 
or equal to k and A their corresponding lower weights obtained in the 
usual way by subtracting multiples of simple roots en from A until 
the lowest weight is reached. The specific choice of the groups ensures 



that these theories possess a mass gap [16]. Of particular interest is the 
dimension of the perturbing operator, that is the second term in (|l|), 
which is associated to the conformal dimension related to the adjoint 
representation A = ip = the longest root of G . We obtained the unique 
dimension A(ip, 0) = hj (k+h). Since A(vp, 0) < 1 for all allowed values 
of k, this perturbation is always relevant in the sense of renormalisation. 
To compute the conformal dimensions A(A, A) of the entire primary 
field content is in general quite a formidable task. We report here only 
the SU(3)2/U(1) 2 example with c = 6/5, which we treat in detail in 
several different approaches. We obtained j| 12 fields with dimension 
1/10, 6 fields of dimension 1/2, and one field, i.e. the perturbing 
operator of dimension A(^,0) = 3/5. 

Having established that the theory defined by (||) constitutes an 
integrable massive quantum field theory suggests the existence of a 
factorizable scattering matrices which can be associated to it. Indeed, 
based on the assumption that the semi-classical spectrum is exact, the 



S-matrix elements have been conjectured in [17] and verified by means 
of the bootstrap principle for HSG-models related to simply laced Lie 
algebras. Instead of just stating these results, I would like to provide a 
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slightly broader context and present a recently found type of scattering 
matrix of which the one of [jl^] is a special case. 

3. Colour valued S-matrices 

For this purpose I want to recall briefly the key features of step i) of 
the bootstrap approach. The fundamental observation, on which all 
further analysis relies upon, is that integrability, that means here the 
existence, one does not need to know its explicit form, of at least one 
non-trivial conserved charge, in 1+1 space-time dimensions implies the 
factorization of the n-particle scattering matrix into a product of two 
particle scattering matrices 

z» n (e n ) . . . z^ie x ) |o) out = J] s^ie^zM . . . z,je n ) |o> in . 

l<i<j<n 

(3) 

It is usually convenient to parameterize the two-momentum p by the 
rapidity 6 as p = m(cosh 9, sinh 9). We abbreviate 9 ab '■= &A — 0B- The 
Z^{9) are creation operators for stable particles of type fx with rapidity 
9, which obey the Zamolodchikov algebra. The basic assumption of the 
bootstrap program is now that every solution to the unitary-analyticity, 
crossing and fusing bootstrap equations^ 

Sab(O) = SBAi-9)- 1 = S BA (in - 9), ]J S Dl (9 + m m ) = 1 , 

l=A,B,C 

(rji € Q are the fusing angles encoding the mass spectrum and the anti- 
particle of A is A), which admits a consistent explanation of all poles 
inside the physical sheet (that is < Im# < ir), leads to a local quan- 
tum field theory. There exists no rigorous proof for this assumption, 
however, it is supported by numerous explicitly constructed examples. 
We may now exploit this principle and construct new solutions (2| . For 
this purpose we supply each particle type with an additional substruc- 
ture and identify each particle by two quantum numbers, i.e. A = (a, i), 
such that the scattering matrices are of the form S l J b (9). We associate 
the main quantum numbers a, b to the vertices of the Dynkin diagram of 
a simply laced Lie algebra g of rank t and the colour quantum numbers 
i,j to the vertices of the Dynkin diagram of a simply laced Lie algebra 
g of rank i. Noting now that for many (not always though) theories 
the scattering matrix factorizes into a so-called minimal part S^^(9), 

1 For the purpose of this talk I suppose that there is no backscattering in the 
theory such that the Yang-Baxter equation constitutes no constraint. 
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which satisfies §) by itself, and into a CDD-factor S^ D {6 , B (p)) 
which depends on the effective coupling constant B{(3) containing only 
unphysical poles, S A b(0) = Sfg (0) S^ D (6 , B (/?)), we may include the 
colour structure by defining the interaction as 

Allowing now g,g to be simply laced algebras, we found [Q] as solution 

= e^^W+J^ f (2coshf -7) i .(2co S hf-j) a - 1 e-«^)_ (g) 

The theory possess various free parameters, the £ mass scales of the 
stable particles and the I mass scales of the unstable particles char- 
acterized by the resonance parameters cxij = (1 < i,j < £) with 
the additional constraint that i and j are connected on the g-Dynkin 
diagram. 

As special cases of these models we have the g| Ai-theories and the 
A n |g-theories, which coincide with the minimal affine Toda theories 
and the g n +i-HSG-models, respectively. S-matrices which allow the 
"colour algebra" g also to be non-simply laced were recently proposed 
in [18]. The generalization which also admits g to be non-simply laced 



still remains an interesting open problem. In that case the outlined 
factorization into minimal and CDD-factor does not take place anymore 
and one therefore has to abandon the described construction principle. 

As already mentioned as particular example, which we want analyze 
in more detail, we consider the A1JA2 = S'[/(3)2-HSG-model with 

S±± = -l , S± T (6>) = ±tanh(0±CT-i7r/2)/2. (7) 

Since there is no pole in the physical sheet, fusing processes do not 
take place. An important point to note here, to which we frequently 
appeal, is that for the limit a — > 00 the S-matrices S±^ tend to one 
and the whole theory decouples into a direct product of two thermally 
perturbed Ising models. 



4. Description of unstable particles 

Since the states in (||) are asymptotic, it is clear that unstable particles 
have to be described by different means, even when they possess a very 
long life time. As is familiar from standard quantum mechanics, one 
may introduce a decay width by complexifying the mass of a particle. 
As is for instance explained in the same prescription can be taken 



Kiev2000.tex; 1/02/2008; 12:01; p. 5 



6 



Andreas Fring 



over to quantum field theory. The description of an unstable particle of 
type c may be thought of as adding a decay width to the physical 
mass of a stable particle, such that the S-matrix as a function of the 
Mandelstam s-variable has a pole at s = = (Mg — iTc/2) 2 . Trans- 
forming from the s to the rapidity plane and describing the scattering 
of two stable particles of type a and b with masses m a and m& by S a b(0), 
the resonance pole is situated at Or = a — ia. Identifying the real and 
imaginary parts of the pole then yields 

M d — = m a + m b + 2m a mb cosh a cos a , (8) 

McTc = 2m a rrn ) sinh |er| sina . (9) 

Whenever Mg 3> Tg, the quantity Mg admits a clear cut interpretation 
as the physical mass and the relations (ph and fflf) acquire the form 



which is usually referred to as Breit-Wigner formula [20|. Since this 
assumption is only required for interpretational reasons we will not 
demand it in general. Eliminating the decay width from (|8|) and @, 
we can express the mass of the unstable particles Mg in the model as a 
function of the masses of the stable particles m a , nib and the resonance 
parameter a. Assuming a to be large we obtain 

M| ~ m a m b {l + cosct)/2 e H . (10) 

The occurrence of the variable me'' 7 '/ 2 is interesting, since it was 
introduced originally in |2l]] in order to describe massless particles, i.e. 
one may perform safely the simultaneous limit m — > 0, a — > oo, and one 
might therefore be tempted to describe flows related to (10) as so-called 



massless flows. In fact, in |I| we introduced the variable explicitly as a 
formal parameter. 

One should note, however, a few essential differences between the 
roaming or staircase models || and the HSG models. First of all the 
decay width of the unstable particles enters in a different manner. 
Whereas in the HSG-model they are introduced by a rapidity shift 
6 — > 9 ± a of a parity broken theory, in the staircase models they enter 
through an analytic continuation at the self-dual point of the effective 
coupling constant B, i.e. B — ► 1 ± ia/ir, of scattering matrices associ- 
ated to affine Toda field theories. Concerning the analytic continuation, 
the point B = 1 is quite special, since for this point the combination 
B(B — 2) remains real. The other distinction between the two classes of 
models is the origin of the staircase pattern observed in the scaling func- 
tions of the models (see below). For the HSG-models one may associate 
the steps directly to the energy scale of the unstable particles, which 
is not possible for the staircase models. Furthermore, the HSG-models 
admit a well-defined Lagrangian. The condition Tg <C Mg corresponds 
to the semi-classical limit. 
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5. Thermodynamic Bethe Ansatz 

Originally the TBA was formulated in the context of the non-relativistic 



Bose gas by Yang and Yang [22] in order to extract various ther- 



modynamic quantities. Thereafter it was extended [|23J to relativistic 
quantum field theories whose scattering matrices factorize into two- 
particle ones. In this context the approach allows the construction 
of a scaling function which reproduces some characteristic features, 
to be specified in more detail below, of the model and in particular 
various quantities of the underlying ultraviolet conformal field theory, 
in particular the effective Virasoro central charge and sometimes also 
the dimension of the perturbing operator. 

It is well-known that in a relativistically invariant theory one can 
in general not use a wavefunction formalism, because of the creation 
of real and virtual particles. However, in configuration space one may 
assume to have regions in which one regards the particles as free and 
describe them formally by a function ip(xi, . . . , xjv), N being the total 
number of particles. Two adjacent free regions are then connected by an 
S-matrix. Taking now the particle A on a trip around the world, a circle 
in this case, the formal wavefunction of A picks up the corresponding 
S-matrix element as a phase factor when meeting another particle. For 
the case at hand we have to distinguish, whether the particle is moved 
clockwise or counter-clockwise along the world line, due to the breaking 
of parity invariance. Proceeding this way, we obtain two types of Bethe 
Ansatz equations 

e iLM A su^e A "Q Sab{6ab) = e -iLM A sinh dA JJ fi^^) = lj (H) 
B^A B^A 

with L denoting the length of the compactified space direction. These 
two sets of equations are of course not entirely independent and may 
be obtained from each other by complex conjugation and noting that S 



is Hermitian-analytic [ 24 1 . Unlike as for instance in the talk of Slavnov 
(see this proceeding) , we do not use the equations (|r[|) to find allowed 
momenta, but instead we take the logarithm of these equations and 
carry out the thermodynamic limit in the usual fashion [23], i.e. letting 
the size of the system go to infinity L — > oo, while keeping the ratio 
L/N finite. This way we obtain the following sets of coupled non- linear 
integral equations 

Y, B §AB*L%(6,r) = rM A cosh + ln(l - exp(-L±(0, r))) . (12) 

Here the rapidity convolution of two functions is denoted by f*g(8) := 
/ d9'/2ir f(9-9')g(9'). We also re-defined the masses by M A -> M A fm\ 
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keeping the same notation. The parameter r = miT -1 is the inverse 
temperature times the overall mass scale of the lightest particle. The 
kernels in the integrals carry the information of the dynamical and 
statistical interaction of the system and are given by 

$AB(d) = <5>ba{-9) = _i dlnS d A ^ _ 2ttW(0). (13) 

For generic values of qab the statistical interaction is of Haldane type 
[p5fl . The special cases gAB = and g^B = &AB correspond to bosonic 
and fermionic statistics, respectively. Recently also a formulation for 
Gentile statistics has been proposed |2(|. Here we will restrict our atten- 
tion to the choice gAB = $AB- As very common in these considerations, 
we could also parameterize the function LAO) = ln(l + e~ e AW) by the 
so-called pseudo-energies e A (9). Then it follows from the properties 
of the TBA equations that (-\{0) = e^(— 0). In comparison with the 
parity invariant case the main difference is that we have lost the usual 
symmetry of the pseudo-energies as a function of the rapidities, since 
we have e^(0) ^ tjjfi) m ^ ms case - The main task of the TBA is now 
to solve the equations (|l2^) for the L's. Thereafter, one may compute 
the scaling function 

c(r) = ^Y, MA J d9 cosh 9 (L^{6) + L+ {6)), (14) 
whose key property is 

lim c(r) = c eff = c - 24A' , (15) 

r—*0 

with c being the Virasoro central charge and A' the lowest occurring 
conformal dimension of the underlying ultraviolet conformal field the- 
ory. We recall that for unitary models the values of c and c c fr coincide. 



Due to the non- linear nature of the equations (12), there exist hardly 
any exact analytical solutions for the full TBA-equations, albeit it is 
possible to find various analytical approximations, e.g. p7|]. When one 



is only interested in c fj, it is possible to solve (12) and (0) fully an- 
alytically. Noting that the TBA-equations for the g|g-theories develop 
constant regions in this limit one may write (12) in this case in the very 
symmetric form 

UiQlf ab + U(Qif' = (Q l a) 2 ■ (16) 

b=l j=l 

Here we introduced the variable Q l a = Yib^i-^l)^^ ^ or nm t e 
values of the resonance parameter the effective central charge is then 
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expressible as 



-eff 



'h 



6 



h + h 7T 



a=l i=l 



CIO 



(17) 



with £(cc) = YlnLi x n /n 2 +lnxln(l— x)/2 denoting Rogers dilogarithm. 
Remarkably, besides the numerical solutions, for many combinations of 
g and g, the relations ([If]) have been solved analytically [gg, ^9|, [l], |3(| 
in such a way that the Q^s acquire the form of certain Weyl-characters. 
Recalling that the A n | g-theories= gn+i-HSG-models we note that the 
first part of(||) and ([17]) obviously coincide. 

From a purely mathematical point of view the system of equations 
(16) and second part of (|17]) is an exceptional situation whenever c er r 
is rational, in which case it is referred to as accessible dilogarithms 
(for a review see e.g. [ j3l| l and references therein). Remarkably, the 
same system of equations may be obtained by employing techniques 



originally pursuit in [32] 




su(3), 


m 1 /m 2 =' 


,(T=0 






,o-=65 


su(3) 3 


m 1 /m 2 =' 


,cr= 50 



log(rfi) 



(b) 



Figure 1: Scaling functions for su(3)^ , k — 2, 3, 4 from the thermodynamic Bethe ansatz as a function of 
the variable logr/2 with r — m\T 1 at different values of the resonance shift and mass ratios (a), (b). 
Scaling function from the correlation function {©(r)©(0)) for su(3)2 (c) and RG-fiow of the conformal 
dimension of the local operator Cq'q (d) as a function of 2 log r. 



Kiev2000.tex; 1/02/2008; 12:01; p. 9 



10 



Andreas Fring 



In [|J we analyzed the full TBA-equations for the Afc_i|A2-theory= 
SU(3)fc-HSG. The scaling function may be approximated in this case 

where is the smallest mass of an unstable bound state which may 
be formed in the process (a, i) + (6, j) — > (c, k). This analytical approxi- 
mation is confirmed by our full numerical solution as depicted in figure 
1(a), (b). 

For the different algebras treated, we note that the different plateaus 
in the scaling functions are characterized by the mass scales of the un- 
stable particles as described in ( |i~8| ) . By choosing the masses of the two 
particles to be very different we also confirm in part (a) the decoupling 
of the theory into a direct product of two thermally perturbed Ising 
models. Exploiting the observation [||, ||, |(| that the plateaus are also 
realizable as cosets, it is also possible to find analytical descriptions 
similar in spirit to (|l6|) and (|P7|) for them ||30||. 



6. Form Factors 

To use form factors as auxiliary objects in order to compute correlation 
functions is an approach originally pioneered by the Berlin group in the 
late seventies Form factors are tensor valued functions, represent- 
ing matrix elements of some local operator 0(x) at the origin between 
a multiparticle in-state and the vacuum, which we denote by 

F^-^(et, . . . A) := <0|O(0)|Z Ml (0i)Z M2 (0 2 ) . . . Z, n (9 n )) in . (19) 

Here the Z At (0)'s are the same type of creation operators as in Eq. 

One way to determine the form factors, that is step ii) in the boot- 
strap program, is similar in spirit to the computation of the S-matrix, 
namely one considers the consequences various physical concepts have 
on these matrix elements. This way one can set up a system of con- 
straints, which, like in the case of the scattering matrix, turns out to 
be so restrictive that they lead to their explicit determination. We shall 
now only recall these constraints and refer for a more systematic and 
rigorous derivation of them to [34]. 



As a consequence of CPT-invariance or the exchange of two opera- 
tors Z^{6) one obtains 

F n ^ 3 (...6i,9j...)=F n (. . . 9j, 0i . . O-S^At ■ (20) 
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The analytic continuation in the complex #-plane at the threshold cuts 
when = 2ni together with crossing leads to 

f c?K...m„ {9i + 2nij ^ 0n) = F ow..^ (02> _ _ j ^ 9i) (21) 

Since we are describing relativistically invariant theories, we expect for 
an operator O with spin s 

F^-^(e 1 + \...,e n + X) = e sX F^- fi -(e 1 ,...,e n ) , (22) 

with A being an arbitrary real number. For a form factor whose first two 
particles are conjugate to each other we have the so-called kinematical 
pole at in, which leads to a recursive equation relating the (n-2)- and 
the n-particle form factor 

-i _Res F^f^-^ie.+m^e,,...^) = 

(1 - w JJ" =1 WM) F°^-^(e v . . .,e n ), (23) 

with lo being the factor of local commutativity and p, the anti-particle 
of n. We restrict our initial considerations to a model in which stable 
bound states may not be formed and therefore do not need to report 
on the so-called bound state residue equation, which relates an (n-1)- 
and an n-particle form factor. 

The equations (^)-((2^) allow already the determination of several 
solutions. However, since these equations are completely independent 
of 0(x), one needs further arguments to establish which type of oper- 
ators one deals with. This means we are in a quite different position in 
comparison with the computation of matrix elements involving some 
well-known operator and Bethe states as described for instance in the 
talk by Kitanine (see this proceeding). There are however additional 
properties, which do involve the nature of the operator. In general 
they make use of their asymptotic behaviour, like the momentum space 
cluster property. It states that whenever the first, say k, rapidities of 
an n-particle form factor are shifted to infinity, the n-particle form 
factor factorizes into a k and an (n — K)-particle form factor, which are 
possibly related to different types of local operators 

7lVf(0i, • • • ,0 n ) ~ F? . . . ,6 K )F°" K (9 K+1 , . . . ,0 n ) . (24) 

We have introduced here the operator T^ b := lim^oo Y\p=a which 
is composed of the translation operator acting on a function of n 
variables as T a A /(0i, . . . , a , . . . , 9 n ) ^ /(0 1; . . . , a + A, . . . , 9 n ). For 
the purely bosonic case this behaviour can be explained perturbatively 
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by means of Weinberg's power counting theorem, see e.g. [33, 35]. For a 



similar restricted situation, there exists also a short heuristic argument 



which provides some form of intuitive picture of this behaviour [37] by 
appealing to the ultraviolet conformal theory. However, this argument 
is based on various so far unjustified assumptions, which need further 
clarification. In addition, up to now it has never been made use of. It 
is very interesting that when ( f24| ) is read from the r.h.s. to the l.h.s. it 
may been seen a kind of fusion (7^ the fusing mentioned in section 3 
which refers to the particles) of the local operators. Therefore, in order 
to shed more light on step iii) of the bootstrap program, it would be 
important to put the property (p4|) on the same solid ground as the 
equations ©-(H)- 

Having determined the form factors one may carry out step iv) of 
the bootstrap program and expand the two-point correlation function 
as 



00 „ 
(O(r)O'(O)) = E E / 

71=1 



d6\ . . . d6 n _ r E 
n!(2vr) n 6 



x F^-^(e 1 ,...,e n ) (F°'^-^(9 l ,...,e n ) ) (25) 



where E = X/ILi m fj.i cosh#j is the sum of the on-shell energies. Pro- 
vided we can compute (]25| ) we may use the two-point correlation func- 
tion in order to determine the RG-flow of the Virasoro central charge 
c and the conformal dimension of a primary field. According to the 
c-theorem of Zamolodchikov we have 

c{r) = - dss 3 (9(s)G(0)) (26) 
2 ./,• 



oo 



71=1 /Jl...M n 



d0\ . . . d9 n _ rE 
n!(2vr) n 6 

2 (6 + 6rE + 3r 2 E 2 + r 3 E 3 ) 
2E~ A 



(27) 



where is the trace of the energy momentum tensor. The flow of a 
conformal weight of a primary field was proposed in || 

Aw =-2(c5(o))r <i5s<e(s)o(o)> <2s) 

no ° d9 l ...d6 n {l + rE)e- rE 



n!(2vr)" 2 (O(0)) E 2 



-E E 

n=lMl-Mn 

xi^-^O?!, ...,9 n ) (if 1 ^---^!, ...,0 n )Y (29) 
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as a straightforward generalization of the A-sum rule [37]. The con- 
formal dimension of a primary field may also be obtained by a di- 
rect analysis of the well-known ultraviolet behaviour of the correlation 
function 

lim (O(r)O(O)) ~ r~ 4A ° . (30) 

r— >0 



The main difference between the two possibilities is that ( |30| ) leads 
to far less conclusive values as (|29|). Whereas ( |30[ ) requires in general 
an a priori knowledge of A , and serves mainly to confirm this, the 
expression ( |29| ) can just be evaluated by brute force. 

I will now present the explicit solutions found for equations (|20|)- 
© for the SCf(3) 2 -HSG model |, |, |]. A step towards a more group 



theoretical understanding was done in |6[, where these results were 
generalized to the SU(N)2 case. 

Labeling an operator by four quantum numbers fi, v, r, r 1 we proved 
in || that the general n-particle solution reads 

O^AM+M- 0^" v ,\M+M- 
rp t,t" In a \ tt t.t" i , 

r 2s+T,2t+r' {"If ■ ■ ■ i °n) ~ IJ 2s+T,2t+T l ufc!l " /i 2s+T,2t+r' 

i<j 

W.l.g. we assumed here a particular ordering by starting with 2s + r 
particles of the type fi = + followed by 2s + r' particles of the type 
fj, = — , collected in the sets M = {±,...,±}. Once these expres- 
sions are known, all other form factors related to it by permutations of 
the particles may be constructed trivially by exploiting (p0|)-(21). The 
functions for all combinations of the /x's are 

F±±(#) = -i/2tanh^exp(^/2) (32) 

m(l T l)±6> _G dt sm 2 ((in-e^a)±-) 
p±T(Q~j — 2Ie 4 tt JO t sinhtcoshi/2 (33) 

with G being the Catalan constant. The (t + s) x (t + s)-matrix 
/ \ f &tr ■ -\ , , 1 < i < i 



has as its entries elementary symmetric polynomials depending on dif- 
ferent sets of variables. We use the notation when they depend on 
the variable x = exp 6 associated to the sets and a to indicate that 
all variables are multiplied by a factor ie~ a . The overall constant was 
computed to 

ZT°r'y _ ■s(2t+t' '+u+2)r > s(2s-2t-r' -1+2-r) e sa(2t+r')/2 tt°t',t> /or^ 

n 2s+r,2t+T> — % Z e n r,2t+r' \ 60 ) 
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0"'", 

where the value of H T 2t+r' 1S fixed by the lowest non-vanishing form 
factor. To determine this constant requires additional arguments. That 
these expressions are indeed solutions of (|20|)-(23) was rigorously proven 



in ||. Alternatively we showed that the following cluster properties hold 

rf\ P M,0 p 0,0 p A*+£(l-2jti),0 ,„„x 

J l,2K+£<l r 2s+T,2t+T' - r 2/t+5,0- r 2s+r-2K-5,2i+r" 

I l,2K+£<l r 2s+r,2t+v r 2K+£,0 r 2s+T-2K-£,2t+v> V°'J 

rf\ r 0,u p 0,H-£(l--2i/) p 0,0 

7 n+l-2K-5<m,n- r 2s+T,2t+r' - r 2s+T,2t+-r'-2K-£- r i 2K+£' \°°J 

q—X F n,v j?n,v ^0,0 

' i n+l-2K-g<m,n- t 2s+/x,2t+T' J 2s+/i,2i+"r'-2K-£ J - 0,2k+? ' V 03 / 



where we simplified the notation in an obvious way. The equations (36) 



(39) may be seen on one hand as a consistency check, since we recover 
the solutions constructed directly, and also as a construction principle, 
because from one particular solution we may derive various other ones, 
i.e. we simulate the "fusion" process. 

Since the form factors of the trace of the energy-momentum tensor 
occur in various applications it is instructive to use it as an example 
to give the idea how more explicit formulae look like when all the 
quantities are assembled together 

F 2s,2t = <ri(xi, . . . ,x n )cTi(2;j" 1 , . . . ,x~ 1 )F 2s 2 f t . (40) 

To fix the initial condition is in general also still a difficult issue, but we 
may solve this problem here easily by appealing to the fact that when 
the n-particle form factors involve only one type of particle they have 
to coincide with the ones for the thermally perturbed Ising model. The 
only non-vanishing form factor for the trace of the energy-momentum 
tensor is 

Ff{0) = -27rim 2 sinh(0/2) . (41) 
Explicit expressions for 4,6-particle form factors are for instance 

F e|++— -7rml e (03i+fl42)/2(2+S. < .cosh(e <j )) yr M j(ft , , 

r i — 2cosh(6»i 2 /2)cosh(6»34/2) H r min K a *j) 

i<j 



F e\ +++ +- _ 7rm 2 (3 + Ei<j cosh(%)) 



Kj 
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The matrix A is in this case of dimension (t + s — 2) x (t + s 
reads explicitly 

0\ 




2) and 



A 



e 



0"r 







0"l 



°3 



(Jo 



V 



















£7 



2s-l 






(44) 



(-l)^2 t -l / 



Exploiting the integral representation for the symmetric polynomials 
<T fc (xi, ...,x n ) = t^i § U \ z „ dz k+1 l\k=i(. z + x k)i it is easy to express 



the determinant of A® as a multidimensional integral representation 
det^ e = (-l)( s+1 )* 4) dm ... i </»,_, jj </r ... I dr„ 



t-1 2s 
2-2s-2j 



3=1 



Y[(uj + Xi 



1 2s+2t 
2-2*-2j 



■Js-1 



i=l 



n 

3=1 



>< n k 2 - u i) n ^ 

l<i<j<t-l l<i<j<s-l 

where we abbreviated <f> dz = (2m) t 



n < 

i=l+2s 

8-1 t-l 



Vj + Xi 



(45) 



j=l 1=1 

dz with q being a positive 



real number. These type of integrals differ from the ones provided in 
[|34|1, where the integration is carried out in the 6 rather than the x- 
variables. To manifest the precise relation between the two formulations 
remains an open issue at present. 

Indicating by a superscript the number of the highest n-particle form 
factors taken into account in ( |27| ) and (p9|), we computed by means of 
a Monte Carlo integration 



c(0) (2) = 1, c(0) (4) = 1.197..., c(0) (6) = 1.199... , a < oo (46) 
lim c(0) = 1 . (47) 

IT— >CO 

The latter equation is exact, since we have linio-^oo i^+^f ~ e ~(. t + s ~ 1 )' 7 ^ 

from which we deduce that only .F®' ' 2 and i 7 ®' 2,0 are non- vanishing, 
which allows an analytic evaluation. We obtained [Q] 

A°o,o(o)( 6 ) = 0.1004 cr<oo. (48) 
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The entire RG-flow of c and A -° is depicted in figure 1. We observe 
two crucial features: First of all, by comparing part (a) and (c) we see 
that the c-flow determined from is qualitatively identical to the 
one obtained from the TBA. Secondly by comparing part (c) and (d) we 
observe that whenever the S'C/(3)2-HSG model decouples into a direct 
product of two thermally perturbed Ising models, the o ' o -field decou- 
ples into two fields with conformal dimension 1/16, i.e. the disorder 
operator of the Ising model. We confirm these findings alternatively by 
a direct analysis of ( |30|) . A similar computation y] was carried out also 
for several other solution (31). Unfortunately, we did not find solutions 
corresponding to all primary fields, in particular we did not find any 
solution corresponding to A = 1/2. 



7. Conclusions and open issues 

The main conclusion which can be drawn from our analysis concern- 
ing the specific status of the HSG-models (g|g-theories) is that the 
scattering matrix originally proposed in [jl7| is certainly consistent and 
can be associated to a perturbed gauged WZNW-coset. This is based 
on the fact that we reproduce all the predicted features of this pic- 
ture, namely the expected ultraviolet Virasoro central charge, various 
conformal dimensions of local operators and the characteristics of the 
unstable particle spectrum by means of the TBA as well as the form 
factor analysis. 

The key statement concerning a comparison of the two different 
methods is that they yield the same qualitative behaviour of the scaling 
function. With regard to the efficiency of the methods one should point 
out that in the TBA approach the number of coupled non-linear integral 
equations to be solved increases with the number of different particles, 
which means the system becomes extremely complex and cumbersome 
to solve even numerically. Computing the scaling function with the help 
of form factors only adds more terms to each n-particle contribution, 
but is technically not more involved. The price we pay in this setting is, 
however, a slower convergence of (p5|). With regard to other quantities, 
that is the flow of the conformal dimensions, one should note that in the 
TBA-approach it is not known how to solve this problem and even at 
the fixed point itself one is, so far, only able to determine the conformal 
dimension of the perturbing operator with considerable effort for some 
particular theories. In the context of form factors there is no obstacle 
to the computation, once the solution is known. 

Concerning the status of the bootstrap program, one should say 
that a very rich structure can be associated to step i). In particular, 
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there exist closed general group theoretical expressions, which allow 
a relatively universal understanding as opposed to model dependent 
statements. Nonetheless, a complete classification is still outstanding. 
The status of step ii) is still less satisfactory since computable closed 
formulae exist only in few cases and one is still far from a closed group 
theoretical formulation. Step hi) also still requires some considerations, 
since up to now it is not even possible to identify properly the operators 
related to primary fields (we did not find various dimensions and can 
not undo the degeneracy), not to mention all the remaining operators. 
Further study of the cluster property will certainly provide advance in 
this direction. The computation of correlation functions, i.e. step iv) can 
in principle be carried out numerically, but of course a better analytical 
understanding is highly desirable. Our results || also indicate that 
the "folkloristic belief" of the fast convergence of the series expansion 
of (^5|) has to be challenged. In fact, for large values of N, this is 
not true anymore. It would be highly desirable to have more concrete 
quantitative criteria at hand. Needless to say that it would be nice to 
increase amount of technical tools available to carry out v) and enrich 
the spectrum of non-perturbative arguments. 
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